Colloidal particles at a nematic-isotropic interface: effects of confinement 
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When captured by a flat nematic-isotropic interface, colloidal particles can be dragged by it. As 
a result spatially periodic structures may appear, with the period depending on a particle mass, 
size, and interface velocity p]. If liquid crystal is sandwiched between two substrates, the interface 
takes a wedge-like shape, accommodating the interface-substrate contact angle and minimizing the 
director distortions on its nematic side. Correspondingly, particles move along complex trajectories: 
they are first captured by the interface and then 'glide' towards its vertex point. Our experiments 
quantify this scenario, and numerical minimization of the Landau-de Gennes free energy allow for 
a qualitative description of the interfacial structure and the drag force. 



I. INTRODUCTION 

Precise manipulation of tiny particles in liquids is 
a rapidly developing direction of a modern technol- 
ogy. Manufacturing of e-papers and electrophoretic dis- 
plays 13, , separation of bacterial species and living 
cells J 4 , trapping of DNA and polymerparti- 

cles 0, H, 05 growth of photonic crystals [To, lU, are 
only a few examples revealing its scientific and techno- 
logical importance. 

In order to move, organize, or separate colloidal parti- 
cles several techniques have been suggested, among them 
the drag of the micro-particles by a moving nematic- 
isotropic interface 1]. Briefly, due to the interfacial and 
surface tensions as well as long-range distortions of the 
director field, particles can be captured and subsequently 
dragged by the moving interface. Matching the speed of 
the interface, particle size, and elastic properties of the 
liquid crystal, one can move particles of specified radius 
and control their spatial distribution in the cell. 

In our previous work 1] we consider the situation when 
the interface between the nematic and isotropic phases is 
flat. In a slab geometry, however, the interface bends 
accommodating the contact angle between the nematic 
phase, the isotropic phase, and the substrates, and min- 
imizing the director distortions on the nematic side of 
the interface 12, 13]. As a result, the meniscus of the 
isotropic phase extends into the nematic phase; this of 
course influences particle trajectories in the interfacial 
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FIG. 1: Experimental setup. 



region. In this paper we study the behavior of colloidal 
particles captured by a wedge-shaped interface. We first 
examine the shape of the interface and the director dis- 
tribution in its vicinity. Then we analyze how these fac- 
tors influence the trajectory of a particle in the inter- 
facial region. We conclude that there is an additional 
force on a particle acting towards the vertex of the in- 
terface. Finally, we compare our experimental results to 
the estimates based on the minimization of the Landau- 
de Gennes free energy. 
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FIG. 2: Typical optical images of the extended interface. The 
microscope objective lens is focused on a) bottom substrate; 
b) vertex line, cell bulk; c) top substrate. 



FIG. 3: Snapshot of the moving interface. Hot stage cooling 
speed 5deg/min. Orientations of the polarizer (P) and ana- 
lyzer (A) are marked by the white arrows. Rubbing direction 
is parallel to the horizontal side of the pictures. 



II. EXPERIMENT 

Spherical polystyrene particles of diaraeter d = 16 jam 
are dispersed in the isotropic phase of liquid crystal 
5CB using ultrasonic shaker. The density of the par- 
ticles, 1.04 g/cm^, is slightly higher than the density of 
pure 5CB, 1.007 g/cm^ J^]. The surface of the parti- 
cles provides planar alignment of the liquid crystal direc- 
tor. 40 /im- thick cells are filled with the mixture in the 
isotropic phase and cooled to room temperature. The 
uniform planar alignment is obtained by rubbing layers 
of polyamide PI2555 deposited on the substrates. 

The experimental setup is sketched in Fig. and is de- 
scribed in detail in Ref. [J. Briefly, one end of the cell 
is placed on a hot stage; the other end is attached to 
a metal plate, kept at a room temperature. The tem- 
perature of the hot plate could be changed with the 
speed O.ldeg/min and is monitored with the accuracy 
0.05 deg. Typical temperature gradients over the cell are 
20—30 deg/cm. Heating/cooling of the hot stage with the 
speed 5deg/min results in the interface movement with 
the speed v ^ 13 //m/s. The system is observed through 
the polarizing microscope. 

Typical optical images of the cell are shown in Fig. [3 
Two homogeneous phases are separated by a rather wide, 
birefringent, region. The width of this region depends 
on the speed of the interface. In a stationary case it is 
about 10 /im, which corresponds to 0.025^ C; at 13 /im/s 
it is about 40 /im. The region occupied by the extended 
interface has a distinct boundary which separates it from 
the nematic side of the cell. We denote this boundary as 
a vertex line. Focusing the microscope lens on the cell 
substrates, we observed two additional lines positioned 
on the right side of the vertex line, in the isotropic phase, 
see Figs. Efa) and (c), contact lines. The vertex line is 
also visible here, however it is less sharp compared to 
the situation when the focus is on the bulk of the cell, 

Fig.ETb). 

These observations allow us to reconstruct the shape 
of the interface. Since the nematic phase of 5CB wets 
polyamide-treated surfaces 0, , the interface bends 



in order to accommodate the zero contact angle at the 
cell substrates [l^, ll^ ll^ • This bending influences the 
director distribution in the cell. The director is anchored 
at the interface with the tilt angle ^ 26.5 deg and 
rather strong anchoring, W ~ 10~^ erg/cm^ 18]. The 
substrates provide planar orientation of the director. To 
match the boundary conditions at the interface and the 
substrates, the director orientation must change from the 
uniform planar orientation far from the interface to the 
tilted orientation at the interface. 

Rotating the cell between crossed polarizers we have 
also noticed that the director is parallel to the rubbing 
direction everywhere except close to the vertex line, see 
Figin^a). This region looked bright between crossed po- 
larizers on the background of the remaining, dark cell, see 
FiglStb). Rotating the analyzer we were able to darken 
the region next to the vertex line. From these observa- 
tions we concluded that, close to the vertex line, the the 
director twists in the xz plane, parallel to the substrates. 
The value of the twist angle was not well reproducible, 
changing from 45 to 90 degrees. 

Director twist can be understood qualitatively, in 
terms of the splay-bend and twist contributions to the 
elastic energy. Next to the vertex line the director has 
significant splay-bend deformation, which is energeti- 
cally unfavorable. It is known that the director can 
relax this deformation by twisting in the perpendicular 
plane [11 S IH ^ . The balance of the elastic energies 
of the twist and splay-bend deformations defines if twist 
occurs or not. Since the twist elastic constant K22 for 
5CB is about two times smaller than the splay and bend 
elastic constants, Ku and i^as, it is energetically more 
favorable for the director to twist than to bend. In fact, 
such an 'escape' of the director has already been observed 
in a similar systeni, where the nematic is in the contact 
with its own melt [13 or with air 13]. 

Now that the shape of the interface and the director 
distribution in the interfacial region are known, let us 
have a look at the motion of colloidal particles. Upon 
cooling of the hot stage, the area occupied by the ne- 
matic phase grows; the interface moves towards the hot 
stage, meeting the particles dispersed in the isotropic 
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FIG. 5: Schematic representation of the capihary force Fc 
and the drag force exerted on a colloidal particle by the 
moving nematic-isotropic interface. Projection F^ of the cap- 
illary force onto the y axis (perpendicular to the substrates) 
pushes the particle towards the vertex point. 



FIG. 4: Snapshots of a moving nematic-isotropic interface, 
separated by the 2s time intervals, which illustrate particle 
attraction to the vertex line before it is captured by the in- 
terface. Two dotted lines show the position of two aggre- 
gates which are stuck to the cell surface and serve as refer- 
ence points. The interface moves slowly from the top to the 
bottom; at the same time, a cluster of two colloidal particles 
(encircled by a solid line) moves towards the interface, from 
the bottom to the top. 



phase. Fast moving interfaces with the velocity above 
18 /im/s are not able to capture the particles. At smaller 
velocities, dispersed in the isotropic phase particles are 
attracted to the interface as soon as they reach one of 
the contact lines, as illustrated in Fig. ^ Immediately 
after that they speed up towards the vertex line, in the 
direction opposite to the interface propagation. Upon 
reaching the vertex line, particles reverse the direction of 
motion and start following the interface. In fact, even 
when the particles are deep in the isotropic phase, there 
is a small force acting on them in the direction opposite 
to the temperature gradient, i.e., towards the interface. 
This force is due to increase of the free energy of a para- 
nematic phase at the surface of the particle when it moves 
into the region with a higher temperature. 

Upon heating, the interface moves away from the hot 
stage, pushing the particles from the isotropic to the ne- 
matic phase. Particles start to interact with the interface 
only at a distance comparable to their own diameter and 
much smaller speeds of the interface, below 4 /im/s, are 
required to capture them. 

To understand particle dynamics we first note that the 
density of the particles is slightly bigger than the density 
of the liquid crystal. As a result, they gradually sediment 
at the bottom substrate, where they meet the moving in- 
terface. If a particle is trapped by the interface, the inter- 
facial tension results in a force on the particle, directed 
perpendicular to the interface. 

Since the particle is moving with the interface it ex- 
periences Stokes drag force in the direction opposite to 
the velocity of the interface. If the velocity of the par- 
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FIG. 6: Snapshots of the moving interface. The three parti- 
cles are attracted to the interface and, because of the vertical 
motion, disappear from the focus (near the NI line, snapshots 
4 and 5). In order to see them in the focus again, we need to 
change the focal distance from z = 90iJ.m to z = 60/im (snap- 
shot 6), where 2; is a vertical distance from the objective lens 
to some arbitrary level. Experiments were made with IGfim 
particles in a 200/im thick cell. 



tide does not change, then the sum of all forces on the 
particle is zero. This means that the elastic plus capil- 
lary forces shall compensate the drag force. However, the 
latter has only the x-component, in the direction of the 
interface velocity. At the same time, the capillary force 
is perpendicular to the interface and has a component 
pointing towards the vertex point, as shown in Fig. [51 
Hence, the particle will move until this component of the 
force vanishes, which is at the vertex point. 

In addition, in order to minimize the elastic free en- 
ergy of the system, the particle tries to occupy the place 
with the strongest distortions of the director [2^. This 
results in an additional force guiding the particle towards 



4 



OS 
■D 



T(X) 







- nematic 


isotropic 







T(0): 



FIG. 7: Sketch of the xy cross-section of the system. Sys- 
tem size Lx X Ly. Boundary conditions: planar ahgnment of 
the nematic director at the particle and cell surfaces. Initial 
conditions: flat nematic-isotropic interface at x = 0. 



the vertex line. Together, these two forces lift the parti- 
cle in the cell, in the direction opposite to the interface 
propagation. In fact, the particle lifting can be directly 
observed by changing the focal distance of the objective 
lens, as shown in Fig. |S1 

To support these arguments, in the following we study 
the behavior of the colloid in the vicinity of the con- 
fined nematic-isotropic (NI) interface in the framework 
of Landau-de Gennes free energy. 



III. THEORETICAL MODEL AND DISCUSSION 

A geometry mimicking the experimental setup is shown 
in Fig. [3 A uniform temperature gradient is imposed 
along the x axis with the 'hot' wall at x = Lx/2 and 
the 'cold' wall at x = —Lx/2. The director at the 'cold', 
bottom and top walls is fixed in the x-direction, and the 
absolute value of the order parameter at these walls is 
fixed to the bulk order parameter of the nematic phase 
at two-phase coexistence. The order parameter at the 
'hot' wall is set to zero. A colloidal particle, which we 
take to be a long cylinder of radius i?, with the symmetry 
axis parallel to the z axis, imposes rigid planar anchoring 
boundary conditions at its surface. This results in the 
formation close to the colloid of two topological defects, 
known as boojums ^2^. The defects are aligned in the 
direction of the liquid crystal alignment, i.e., parallel to 
the X axes. 

The syst em is described by the Landau-de Gennes free 
energy 



HQ} 



ifb + fe)dV 



(1) 



where is the bulk free energy density and /e is the 
elastic free energy density. The minimum of the Landau- 
de Gennes functional F{Q} gives the equilibrium value 
of the tensor order parameter Q. Symmetry arguments 
yield the bulk free energy density 25, 26] 



ft = a(T - T*)TrQ2 _ ^XrQ^ + c [TrQ^] ^ 



(2) 



where the positive constants a, 6, c are taken temperature 
independent. It is convenient to scale out the variables 
by defining 



Qij = 6c/bQij, 



(3) 



It will be understood that such scaling has been car- 
ried out, and we shall omit the over-bars in the text 
below. We also introduce a dimensionless temperature 
r by defining 



a(T - T*) = r6V24c. 



(4) 



The bulk nematic phase is stable for r < tni = 1 
with a degree o f orientational order given by Qh = 
3(1 + ^Jl - 8r/9)/4; Qi,{t > 1) = 0. We model the tem- 
perature gradient in the x direction by assuming that r 
changes linearly with the x coordinate 



Tni (1 + OLx) . 



(5) 



Equation Q implies that the NI transition occurs at x = 
0, and that the unconfined interface would form in the 
yz plane. 

We assume strong planar (parallel) anchoring of the 
director with the particle surface. This is valid if the 
anchoring parameter WR/K >> 1, where W is the an- 
choring energy of the particle surface, and holds for rela- 
tively large colloidal particles or big anchoring strengths. 
We also assume that the nematic phase at the particle 
surface is uniaxial with a scalar order parameter Q^ = 1. 

The elastic free energy density can be written as 



fe 



1^ dQij dQij 
2 ^ dxk dxk 



1 dQij 

2 ^ dxi 



dQi. 



dxk 



(6) 



where the constants Li and L2 are related to Frank- 
Oseen elastic constants by Ku = K^s = 9Ql{Li + 



L2/2)/2 and K22 = 9QlLi/2 and is the bulk nematic 
order parameter. 

We use typical values of the material parameters avail- 
able for 5CB m HI: a = 0.044 x lO^J/m^K, b = 
0.816 X 10^J/m^ c = 0.45 x 10^J/m^ Li = 6 x lO-^^j/^^ 
L2 = 2Li, T* = 307K. The nematic-isotropic transition 
temperature for 5CB is Tni = 308.5 K. 



A. Shape of the interface 

We minimize the free energy (QJ numerically, using fi- 
nite elements with adaptive meshes. The area Lx x Ly 
is triangulated and the order tensor Q is set at all ver- 
tices of the mesh and is linearly interpolated within each 
triangle 27, 28, 29]. The free energy is then minimized 
using conjugate gradients method under the constraints 
imposed by the boundary conditions. 

In the experiment, the 40 jam cells are filled with the 
16 /im colloidal particles. To minimize computational ef- 
forts, we reduce the system size to = 3 /im. In fact. 
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FIG. 8: Cross-section of the cell with the nematic-isotropic 
interface. Color coding shows the n^, out of plane, component 
of the director. The director is practically parallel to the 
z axis close to the interface (twist angle ~ 7r/2, Uz ~ 1). 
Lx = 15/im, Ly = 3/im. The temperature gradient parameter 
a = 0.0036/im~^, which corresponds to ~ 54deg/cm. 



minimization of the free energy in a system of this size is 
a nontrivial task: the smallest length scale involved in the 
problem is given by the thickness of the interface, which 
is of the order of the nematic coherence length at coexis- 

1/2 

tence, ^ = (24Lic/6^) ~ 10 nm. This length scale sets 
the mesh size. Thus, about 10^ mesh points would be 
required for a rectangular uniform discretization of the 
cell. By using adaptive meshes we are able to reduce this 
number by more than two orders of magnitude, which 
makes the minimization procedure computationally fea- 
sible. 

We first have a look at the shape of the interface with- 
out colloidal particle. The chosen anisotropy of the elas- 
tic constants, L2/L1 = 2, favors director alignment par- 
allel to the NI interface, which is in mismatch with the 
boundary conditions at the cell boundaries. As we have 
already mentioned in the previous section, this leads to 
formation of the curved interface. The contour plot of 
the riz (out of plane) component of the director and the 
stream-traces of the in-plane {ux and Uy) components are 
shown in Fig. [SI The interface is quasi- flat in the middle 
of the cell and curved at the horizontal surfaces of the 
box. To avoid strong in-plane splay-bend deformations, 
the director twists at the interface, in a qualitative agree- 
ment with the experimental results. The twist angle at 
the interface reaches 7r/2, i. e. the director is parallel 
to the interface and the substrates. Experimentally ob- 
served twist angle is smaller, of the order of 7r/4. The 
difference between the experimentally measured and cal- 
culated values is due to the different anchoring of the 
director at the interface: for 5CB this angle is about 
26 deg, whereas our model (which takes into account only 
second-order terms in the gradient expansion) predicts 
planar interfacial anchoring. The shape of the interface 
is symmetric, with its vertex point placed in the mid- 
dle between the walls, which is in disagreement with the 
experiment, where slightly asymmetric interfacial profile 
is observed. The asymmetry of the experimental inter- 
faces is most probably due to the non-zero temperature 
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FIG. 9: Free energy U{y) per unit length (circles). Colloid 
radius = 0.25 /xm, system size 9x3 /nm. Temperature gra- 
dient a = O.OOl/xm"^ which corresponds to ^ 15deg/cm. In- 
sets illustrate the shape of the interface for different positions 
of the colloidal particle. Upper panel: x = 0; lower panel 
X = — 0.27/im. Uo{y) (shown in the upper panel, diamonds) 
is the free energy of a colloidal particle in a pure nematic 
phase, at r = tni- Kh ^ 4.93 x lO""*^^ J/m for r — tni — 1. 

gradient across the cell (along the y axis). 



B. Forces on a particle in the interfacial region 

Now let us focus on the interaction of a colloidal par- 
ticle with the NI interface. In our previous work [2^ 
we have shown that a single colloidal particle is attracted 
by a flat interface, with the force which is roughly pro- 
portional to the particle radius. This force drags the 
particles along the direction of the temperature gradient 
(x-axis). However, because of the wedge-like shape of the 
interface, one can also expect an additional component of 
the force, along the y axis. To calculate this component, 
we placed a test particle, modeled as a long cylinder of 
radius R = 0.25 jam with the symmetry axis parallel to 
the z axis, next to the interface and evaluated the free 
energy of the system as a function of the ^/-position of the 
particle, at the fixed x-position. To improve the compu- 
tational efficiency, the system size was further decreased 
to 9 X 3 /im. 

The free energy U{x = const^y) = U{y) per length 
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of the colloidal particle is shown in the Fig. |H1 for two 
values of the x-position of the colloid. The temperature 
gradient parameter a = 0.001/im~^, corresponding to 
~ 15deg/cm. In both cases U{y) has a minimum when 
the particle is in the middle of the cell, at = 0. At 
any other point there is a non-zero force Fy = —dJ^/dy^ 
which pulls the particle towards the middle of the cell 
(upwards, if it is placed at the bottom wall). Therefore, 
even when the interface is not moving, those particles 
which are captured by the interface tend to cluster in the 
middle of the cell, at the vertex line. 

Because of the fixed anchoring conditions imposed on 
both the cell and colloid boundaries, one might expect 
that the particle is repelled from the cell walls even with- 
out the presence of a NI interface. To check this, we have 
also calculated the free energy t/o(^) of a colloidal parti- 
cle immersed in a pure nematic phase (r = tat/, a = 0, 
same cell size and boundary conditions as before). Our 
calculations (see the upper panel of Fig. [UJ demonstrate 
that this force is by the order of magnitude smaller, i. e. 
the presence of the curved interface plays a dominant 
role. 



IV. CONCLUSIONS 

To summarize, our experiments demonstrate that the 
curved shape of the interface and the nonuniform direc- 



tor distribution next to it result in a reach dynamics of 
colloidal particles in the interfacial region. We show that 
the particles can be dragged by the interface not only 
horizontally, across the cell, but can also be moved in the 
vertical direction. This effect can be used to build het- 
erogeneous three-dimensional photonic crystals, where 
the impurities (or holes) shall be positioned in designed 
places 34J. 

The drag effect can also be used to design compos- 
ite materials with complex morphologies. The moving 
interface is able to rearrange the nucleation centers of, 
for example, phase-separating polymer network. During 
polymerization, spatial distribution of these centers af- 
fects the final morphology of the network. This opens up 
an opportunity to create composite materials of extraor- 
dinary properties. 
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